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Abstract 



We study the spatial autocorrelation of energy eigenfunctions V'n(q) cor- 
responding to classically chaotic systems in the semiclassical regime. Our 
analysis is based on the Weyl-Wigner formalism for the spectral average 
Ce (q"'" , q~ , i?) of V'n(q''')'0n(q~)) defined as the average over eigenstates 
within an energy window e centered at E. In this framework Ce is the Fourier 
transform in momentum space of the spectral Wigner function W{'x.,E;£). 
Our study reveals the chord structure that Ce inherits from the spectral 
Wigner function showing the interplay between the size of the spectral average 
window, and the spatial separation scale. We discuss under which conditions 
is it possible to define a local system independent regime for C^. In doing so, 
we derive an expression that bridges the existing formulae in the literature 
and find expressions for Ce(q''', q^, E') valid for any separation size |q^ — q~|. 

PACS numbers: 05.45.Mt, 03.65.Sq,73.23.-b,73.23.Ad 
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I. INTRODUCTION AND MOTIVATION 



One of the key issues in the quantum chaos research is the quest for quantum finger- 
prints of the underlying classical dynamics of generic chaotic systems. In the last decades 
most studies of such kind were dedicated to spectral properties [Q. It is now established 
that, in general, systems with a classical chaotic dynamics are characterized by universal 
spectral fluctuations, which is the so called Bohigas conjecture 0]. Much insight about this 
phenomenon was provided by the semi classical approximation 0] . Specifically, it was shown 
that in the semiclassical regime the energy level density autocorrelation function of a 
chaotic system, evaluated at energy separations encompassing several mean level spacings, 
displays similar statistical properties as those arising from ensembles of random matrices [0 . 
Starting from the random matrix side, advances in proving the connection to the spectral 
fluctuations of chaotic systems were also achieved [§,0. Although a full proof of Bohigas' 
conjecture is not yet available, its domain of validity is fairly established. 

Complementary to the universal view, there is another successful contribution of the 
semiclassical approach to the research in quantum chaos. To every spectrum corresponding 
to a Hamiltonian system there are always deviations from the universal behavior, i.e. system 
specific features. For low dimensional chaotic systems the later are usually nicely explained 
by identifying the system shortest classical periodic trajectories. For instance, knowing the 
actions and the stability of all periodic orbits up to a time T of a given chaotic system, 
the semiclassical approximation explains system specific spectral correlations within energy 
windows AE > h/T. 

The search for an universal behavior in chaotic wave functions is more elusive. In this 
case, one faces a different interplay between the quantum and classical scales where the 
system specific features become important. Actually, the investigation of non- universal 
signatures of the classical underlying dynamics in wave functions, such as scars |T0|, is a 
most fascinating subject to which several theoretical studies were devoted jTlHT^. Whereas 
there are still several open questions to be answered, one can naively say that it is possible to 
separate a "short" universal time from a "long" system specific time regime in the description 
of chaotic wave functions. This is, in a broad sense, the subject of the present paper. 

The simplest statistical measure for chaotic wave functions fluctuations is the two-point 
correlation function C{q^, q^) = ('?/^„(q+)'0*(q~))q, where (. . .)q stands for a local average 
in configuration space and ^/'n(q) is the wave function of the n-th energy eigenstate of a 
given Hamiltonian. To the best of our knowledge this autocorrelation function was first 
discussed by Berry |T5|, although admittedly very similar ideas about the characterization 
of chaotic wave function were already known [p!6| , p!7| . Berry assumed a micro canonical prob- 
ability density in the classical phase space for chaotic quantum states and obtained a simple 
analytical expression for C(q"'",q~). In line with the statistical approach, C was also reob- 
tained by means of supersymmetric techniques in weakly disordered systems [|IB|,|TU|. It was 
also numerically verified by eigenfunctions studies of different dynamical systems [pOH 22 



and experimentally observed in eigenmodes of resonating microwave cavities [^B],^. In all 
such studies the agreement was always very good for spatial separations up to a few wave 
lengths, i.e. "short" distances. By correcting C for contributions of classical trajectories 
for large separations |q^ — q~|, a recent study due to Hortikar and Srednicki ||25| improved 
Berry's formula for C{q^, q~). 
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As a further motivation we add that the understanding of wave function correlations has 
some interesting direct apphcation hke, for instance, in the description of the conductance 
fluctuations of quantum dots in the Coulomb blockade regime. p6|- pO|] With this background 
in mind, we take a fresh look into the question of eigenstates autocorrelations using the 
Wigner-Weyl formalism. By doing it we derive a general expression for wave function spatial 
correlations keeping the semiclassical approximation under strict control. We recover as 
limiting cases Berry's correlation function [|T^] as well as the aforementioned Hortikar and 
Srednicki result We furthermore discuss the important role played by different kinds 



of averaging procedures and open a path for the inclusion of scar contributions. 

The structure of the paper is as follows. In Section |T| we introduce the spectral Wigner 
function W, the object on which this study is built due to is simple relation to the autocorre- 
lation function C{q^, q~). There we also discuss aspects of the semiclassical approximation 
for W which are essencial to understand the different limiting results for C{q^, q~). Section 
PTT| contains the main findings of this study. We show that the semiclassical approximation 
for C is easier to obtain starting from a phase space representation, in particular for smooth 
potentials. We derive expressions for C{q^,q^) depending on the averaging procedure for 
essentially any given spatial separation — q~|. In Section |^ we relate our findings with 
previous analytical and numerical results discussing their validity range. We also include 
three appendices. In Appendix ^ we show the equivalence between our expression for C 
with the one derived in Ref. in a certain limiting regime. Appendices |B| and P are 
technical and devoted to the demonstration of some specific formula appearing in the main 
text. 



II. THE SPECTRAL WIGNER FUNCTION AND ITS SEMICLASSICAL 

APPROXIMATION 



The Wigner function ||^ of an individual energy eigenstate \n) is defined by the Weyl 
transformation 



^^"^ ^n(q^X(q ) exp f-^lq ■ P 



(1) 



where the coordinate = q+ — q^, x = (q, p) is a shorthand notation for the phase space 
point with q = (q+ + q~)/2 and d stands for the number of degrees of freedom of the system. 
iFiom Eq. (0) it follows immediately that, upon averaging over the coordinate space (■ ■ ■)q, 
the inverse Weyl transformation of IV„(x) gives the two-point autocorrelation function 



C(q+, q- 



= (^n(q+)C 

'dp (Wn 




exp 



(2) 



Note that our definition of C is not normalized by a factor (|'?/'n(i")P)q^ as standard. 



The essence of Berry's pioneer work [|15| was to assume that the averaged Wigner function 
of a generic chaotic quantum state \n) of an autonomous system is distributed as a Dirac 
5-function over the surface of energy En, i.e. (iy„(x))q oc 5[if(x) — E^], where -ff(x) is the 
system Hamiltonian. For a Hamiltonian of the form 
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H{x) = pV2m + l^(q) , 



(3) 



this so-called microcanonical probability density in a d-dimensional configuration space leads 
to the well known formula 



C(q+, q-) oc 



Jd/2-i p(q)|q+ - q 1/^ 



p{q)\q+ - q \/n 



d/2~l 



(4) 



where p(q) = 2m[E — V (q)] , = (q^+q~)/2, J,y(x) is the Bessel function of order z/. (The 
formula encountered in Ref. [|T5|] differs from Eq.(^ by a constant due to normalization.) 
Since the above relation does not depend on any system specific features, scaling only with 
the local momentum p(q), it directly reveals an universal behavior of chaotic wave functions. 
More recently, Prigodin and collaborators |]18| microscopically obtained the same result for 
disordered systems using the nonlinear a-model, in a regime resembling quantum chaotic 
systems. Impressive numerical tests of Eq. (^) were presented, for instance in Refs. pT| , p2| , 
reporting on the study of wave functions corresponding to high lying energy levels of the 
two-dimensional {d = 2) conformal billiard. 

In spite of the success of Eq. (H), one proviso ought to be made. Both analytical and 
numerical results only corroborate Berry's conjecture if one regards the average (. . .)q in a 
broader sense. The nonlinear sigma model approach averages over an ensemble of different 
impurities configurations. In addition to the average over q = (q+ -|- q~)/2 covering regions 
encompassing several de Broglie wavelengths, the local averages (. . .)q in Ref. had to be 
taken over all directions of (q"*" — q^) for a fixed value of |q+ — q^| to verify Eq. (^). 

Some time ago Berry |[TT|j3^ formulated a more rigorous approach to this subject, which 
was recently further developed by Ozorio de Almeida It has been shown that Berry's 



original conjecture of microcanonical probability concentration [|I^ is semiclassically verified 
if the average runs over Wigner functions of states belonging to an energy window containing 
several levels This construction is best casted in terms of the spectral Wigner 

function, namely 



(5) 



where the {En\ are the system eigenenergies. The energy smoothing function 5e is for 
convenience chosen as 



5e{E-Er, 



e/tx 



{E-E^Y + e^ 



(6) 



in correspondence to an energy window of width e centered at E. Likewise, we introduce 
the smoothed eigenstate autocorrelation function 



C,(q+, q,E) = /\Y.5,{E- E^) ^„(q+)C(q' 



(7) 



where A is the local mean level density, defined as A = {Y,n ^e{E — En))^ ^ with the average 
taken over the energy levels contained by the energy window e centered at E. (When using 
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the semiclassical approximation, for the sake of consistency, A is obtained from the Weyl 
energy level density, A = 1/pw) As follows from Eq. (|l|), the inverse of Wn directly 
gives V'n(q^) V'nlq") rendering 



C,(q+,q-,E) = a| 



dp 



(2TiTiy 



W\ 



q+ + q- 



, p, £ exp 



p ■ iq 



+ 



q 



(8) 



The advantage of using W is that it provides ways for amenable semiclassical approximations 
in different energy smoothing regimes, allowing for the analysis of Ce at any given spatial 
scale separation. Eq. is the starting point of all results derived in this paper. The 
remaining of this section is devoted to the presentation of the limiting approximations to 
W based on the semiclassical picture of chord and centers, postponing to the forthcoming 
section the corresponding analysis of (^^(q^, q~, 



A. The semiclassical spectral Wigner function 



The spectral Wigner function is related to the Weyl propagator f/t(x), i.e. the Weyl 
transform of the propagator (q"*"! exp(— itif/^)|q~), through [ pJ] , p^ 



W{-K,E]e) = ^Re e""*/^ f/t(x) exp (^Et) . 

nh Jo \h / 



(9) 



The semiclassical approximation of W is directly obtained by inserting in the above equation 
the semiclassical expression for the Weyl propagator p3|, namely 



f/r(x) 



det[ll + A^j(x,t)] 



1/2 



exp 



Sj{x,t) 



i/3j 



(10) 



Here the sum is taken over all classical trajectories with the same traversal time t whose 
phase space endpoints are joined by a chord ^j(x) = x+ — x J centered at x. Figure |I| 
illustrates the phase space structure beneath the semiclassical Weyl propagator. The Maslov 
phase associated with the j-th classical trajectory is given by Pj. In Eq. (|TU]), 1 is the identity 
and Aij is the symplectic matrix (or stability matrix). The latter corresponds to the map 



Sxj' = Aij5xj , resulting from the linearization of the dynamics in the neighborhood of x 



The symbol Sj{x,t) stands for the action, also called center action |33], given by 



Sj (x, t) = j) dq^ - p — J dt H[xj (t)] 



(11) 



where H is the Weyl Hamiltonian, i.e. the Weyl symbol of the Hamiltonian operator. It 
is worth recalling that the Weyl Hamiltonian, H(x), only coincides with the classical one, 
if(x), when the latter is of the form given by Eq. (|^). The first integral at the r.h.s. of Eq. 
(P!T|) is the symplectic area enclosed by the circuit taken along the j-th trajectory connecting 
xj to x^ and closed by the chord — ^j (x) (see Fig. |l]). For autonomous systems the second 
integral is simply the product of the energy Ej corresponding to the j-th trajectory, with its 



traversal time tj. The variation of Sj with respect to the independent variables [Q, leads 



to 
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-dSj/dY> 



dSj I dq , and 



Ej = dSj/dt . 



(12) 



At sufficiently short times for each phase space point x there is only one small chord con- 
tributing to the sum in Eq. (|10|). The short trajectory connecting the chord end points 
has a Maslov phase /3o = 0. In distinction, as t is increased, due to bifurcations there is a 
proliferation of different chords to be summed in Eq. (|TD|). 

The starting point for the semiclassical analysis of W is encountered by replacing Eq. 
(0) into Eq. (g) 



->d+l 



ReE 



dt- 



-et/h 



det[ll + A^j(x,t)] 



1/2 



e^^\-[Sj{^,t) + Et + il3j\ 



(13) 



This formula exemplify the general structure of the semiclassical Weyl representation of any 
quantum object as being given in terms of its classical chord structure in phase space. In 
general, this fact is revealed by the use of the stationary phase approximation to obtain the 
dominant contributions for any observable. This premise will guide our analysis of in the 
next section. In the following, we show how the stationary phase method works in the case 
of the spectral Wigner function given by Eq. (|T^. 
The points of stationary phase are the solutions of 



It 



[Sj{^,t) + Et] = . 



(14) 



This equation fixes the traveling time along the j-th trajectory tj{E) at the energy E, for 
which 



S'j(x 



E) 



(15) 



where the action Sj{:K,E) is the symplectic area corresponding to the ffist integral at the 
r.h.s. of Eq. (|TT]), with the momentum |p| fixed by the energy E. In other words, the 
stationary phase condition selects those trajectory segments which belong to a single energy 
shell C. Thus, for such trajectories, all the chords centered at x have their tips on C. If all 
the stationary points tj (E) are isolated, which is generally the case when the chords centered 



at X are sufficiently separated, we can evaluate W by stationary phase |3^, that yields 



W{^,E;e) 



e-"*^/^Aj(x,E) cos 



h 



+ 7i 



(16) 



where the amplitude is explicitly written as 



M^,E) 



dt • 

detju + A<j[x,tj(E 



dE 



-1 



1/2 



(17) 



and we collected the Maslov phases of the classical contributions in 7j. If the energy shell C 
is closed and convex and x lies inside it there will be always contributing chords to Eq. (|K[). 
To keep the presentation simple we shall only consider convex energy shells in this paper. 
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B. The role of the energy average 



The smoothing e parameter plays an essential role in regulating the convergence of the 
semiclassical approximation for the spectral Wigner function W: as e becomes smaller 
longer classical paths start contributing relevantly to the sum in Eq. (0). It is customary 
to define two characteristic semiclassical scales for e p]T| , |32[] . The first one is the outer scale 
eiarge = ^/Tmin, whcrc Tmin IS the period of the shortest periodic orbit, characterizing the 
typical time to flow around the energy shell. The second one is the inner scale Ismail = ^/th, 
where th is the Heisenberg time defined by the mean level spacing A in the considered 
energy window as th = /i/A. 

Particularly strong contributions to W arise when x is taken in the neighborhood of 
caustics. At such singular points the standard stationary phase approximation is bound 
to fail. When the evaluating point x approaches a caustic of the integrand in Eq. (|13|), 
generically two or more stationary phase points tj{E) coalesce and so do the corresponding 
chords Cj(x). Therefore, we shall also often speak of coalescing chords at caustics. The most 
important kind of caustics influencing will be those at the energy shell itself and the ones 
near periodic orbits on C. As we will show, the first ones are associated with short times, 
whereas the others with the long time dynamics. Correspondingly we distinguish two energy 
averaging regimes in W for points x near the energy shell: (a) e ^ ^large when the signatures 
of all the long trajectories are suppressed and (b) the opposite situation when e < eiarge and 
more trajectories do contribute. Here, to keep the approximation under control, e ^ ^smaii 
is required. 

Let us first consider the case when e ^ ^large- Here, only one short classical trajectory, 
with traversal time to fixed by the stationary phase condition, contributes to W [Q. The 
shorter chord in Fig. |^ serves to illustratate this situation. As x ^ C, to approaches the lower 
integration limit in Eq. (|T3D spoiling the stationary phase approximation. This difficulty can 
be circumvented in the following way. Since the action ^(x, t) is always an odd function of t 
and by changing the cutoff function, exp{—et/h), for an even one with respect to t, we write 



Eq. (|T3D as 

W^(x, E-e)c^-- dt exp -[5o(x, t) + Et] \ . (18) 

2TinJ-oo |det[ll + A^o(x,t)]| ^ 

The resulting integrand displays two stationary phase points located at ±to that coalesce at 
t = as X approaches C. The structure of coalescing stationary points at the origin can be 
obtained by expanding the center action up to third order in t |Q 



5(x,t)^-tH(x)-^t=^i7^i, (19) 

where x is the phase space velocity and Ti the Hessian matrix of the Weyl Hamiltonian, 
both taken at the phase space point x. The integral in Eq. (|TB|) can now be evaluated by 
the uniform approximation method by invoking a suitable change of the integration 
variable. Such transformation, z = z{t), is the one that reduces the integrand phase in Eq. 
(p!8|) to the canonical form z^/3 — j^z (see, for instance, Ref. ||36|). Mapping the stationary 
points ±to into the new ones in z, yields to 7 = — [3S'(x, i?)/2]^/'^. Thus, we obtain 
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W{^,E-e) 



2d+i 



e-"l*°l/%(x,E) 



3^o(x, E) 

2n 



nl/6 



Ai 



'3^o(x,i?)y/'' 



2^ 



(20) 



where Ai{y) is the Airy function |^. This result corrects for a small mistake in Eq. (7.20) of 
Ref. It also contains the relation previously discussed for x taken deep inside the energy 
shell. As the evaluation point x recedes from C, the Airy function argument (in modulus) 
grows very fast. Hence, we are entitled to employ the Airy function approximation for large 
negative arguments and retrieve Eq. (|TB|), provided that for the short trajectory 70 = 0. 

As x further approaches the energy shell C, the spectral Wigner function becomes very 
simple (provided e ^ ^large)- In such situation there is an apparent indeterminacy in the 
amplitude of Eq. ( pO] ) since the symplectic area S'(x, E) vanishes and the amplitude Ao(x, E) 
diverges. In this case it is an accurate approximation to represent the short trajectory by 
the short chord ~ '^ox, and hence the stability matrix A4q becomes the identity. Thus 



5o(x,E) 



12 



4 ^ - H(x)]3/^ 
3 (x Hoi) 



(21) 



and 



dtn 



dE 



1/2 



IxTix 



-1/2 



(22) 



Hence, as it was already shown [^ ,p3 



»'(x,£;;£) 



2 AillPW^ 



x-c \h'±no±\'/'^ i(^'x7i:ox)i/3 



(23) 



In this situation, x — > C and e ^ ^largo for the strict semiclassical regime we easily recover 
the microcanonical probability distribution 



W{x,E;e) ^ 5[H(x) -E] 



(24) 



by recalling that limc^o «~^Ai(y/a) = 5(y). This result does not come as a surprise, it is 
just telling us that we washed out most quantum interference effects reaching the classical 
limit while taking e ^ ^large- It is only by narrowing e that one can explore the rich structure 
of the spectral Wigner function and unveil non trivial quantum features. This discussion 
shall be resumed in a deeper level in the following section, but we can already anticipate 
that Eq. ( P^ is remarkably robust. 

Let us discuss now the case where e < eiarge and x is taken close to the energy shell 
C. Here, one also has to account for pairs of coalescing chords in Eq. (^), schematically 
shown in Fig. |^. These are the short chords corresponding to long trajectories orbiting 
between its tips and winding very closely to a periodic orbit. When x ^ C their traversal 
times and actions, S{'k,E), become degenerate with those of the corresponding periodic 
orbit. Thus, Eq. (pO[) has to be corrected by adding the so called scar contributions to 
the spectral Wigner function, first developed by Berry |rT| and latter refined by Ozorio de 
p^j . The latter formula describes the Wigner scars as a peak of extra intensity 



Almeida 



8 



along the periodic orbits on the energy shell decorated by a fringe pattern. As it was shown 
in Ref. |14] such Wigner scars extend deep inside the energy shell where the spectral Wigner 
function is semiclassically given by Eq. (|T6|) for any arbitrary value of e. Indeed, when 
the action of a periodic orbit is Bohr quantized the contributions of trajectory segments 
for chords whose tips lie on the periodic orbit, can be added in phase to Eq. (|16|). Hence 
the Wigner scars have an enhanced pattern of concentric rings of oscillatory amplitude on 
a two-dimensional surface defined by the centers of all the chords with endpoint on the 
periodic orbit |jT^. Only the edge of this surface corresponds to the domain of the Berry's 
scars formula. This shows that the semiclassical spectral Wigner function is in general 
not restricted to the energy shell, and thus to short chord contributions. In other words, 
we say that the old microcanonical conjecture of Voros and Berry has important 

semiclassical corrections. The distinction between contributions of large and short chords 
will appear again in the study of the spectral autocorrelation function in the next section. 



III. CORRELATIONS OF ERGODIC WAVE FUNCTIONS FOR DIFFERENT 

ENERGY AVERAGING REGIMES. 



In this section we derive a general semiclassical formula for Ce(q"'", q~, £") expressed in 
terms of the system classical chord structure. Our analysis is build on the semiclassical 
approximations for the spectral Wigner function W presented in the foregoing section. 

For any given pair of points in position space, q"*" and q~, the integral in Eq. is 
performed over the entire d-dimensional momentum space. In this study we only consider 
q"^ and q~ within classically allowed regions. Hence, the integration momentum space 
intercepts the energy shell and is naturally divided into a domain located in the interior of 
C and another at its exterior. This is illustrated in Fig. ^. For convex energy surfaces, the 
ones considered here, the spectral Wigner function, whose argument is q = (q"*" + q~)/2, 
exponentially vanishes for values of p in the phase space region exterior to C. Hence, in 
general, the main contribution to the integral in Eq. (^) arises from momenta in the interior 
of C and at its immediate neighborhood. Thus, we are allowed to bound the effective 
momentum integration space in Eq. (H) to the classically allowed momenta. We name the 
so defined integration space the "momentum space associated to q = (q"*" + q~)/2". 

In line with Section ||, first we find the general chord structure that inherits from W . 
This is done by inserting into Eq. (|^) the semiclassical expression for the spectral Wigner 
function given by Eq. (|16D, which is valid for any arbitrary e. The resulting integral can be 
casted as 



C,(q+,q-,E)=A5: 



-etj /h 



/+(q+,q-,E)e 



Jr(q+,q-,E)e-^^^' 



with 



If = 2'' ldpA,ix,E)exp\- 5,(x,E)±p-(q+-q- 



(25) 



(26) 



Now we evaluate if by stationary phase. The stationary phase points pj = Pj(q^,q , E) 
are solutions of 
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d_ 
dp 



+ q" 



q 



(27) 



for Ip. We recall that the variation of ^^(x, E) with respect to the independent variables q 
and p are given in Eq. (P^. Hence, the first member of Eq. (^) is exactly — = — (q^^ ~ 
q~). In other words, the stationary phase condition selects those classical j-trajectories in the 
energy shell whose projected chords in the configuration space match the vectors ±(q"'" — q~). 
This geometrical structure is sketched in Fig. ||, where the projected chords are indicated 
by dashed vectors. In the same figure, the classical j-trajectories are those flowing between 
the intersections of the momentum spaces corresponding to q = q~ and q = q^ with the 
energy surface C. The panel (a) of Fig. ^ corresponds to time reversal symmetric flows, 
whereas (b) represents the cases when this symmetry is absent. The locations of the many 
possible stationary phase points Pj in the momentum space associated to q = (q+ + q~)/2 
are indicated by dots and astericks. While for q^ —>■ q"*" the chords centered at "astericks" 
coalesce to a zero length chord at the energy shell C, the ones centered at "dots" are typically 
large. This allow us to name the chords centered at "astericks" as "short" chords and that 
ones centered at "dots" as "long" chords. 

Provided that the stationary points Pj are sufficiently far apart from each other we can 
safely evaluate Ij' and I J' by the stationary phase method. The latter requires a symmetrized 
Legendre transformation in the phases of Eq. (p6D , which is conveniently expressed by the 
standard textbook action S |Q, with variables in the configuration space, namely 

S,(q±, q^, E) = 5,(x„ E) ± p, • (q± - q^) , (28) 

where = (q, Pj) with q = (q+ + q^)/2 and pj = Pj(q+, q^, i?). We obtain for the non 
oscillatory factor of both and 



2%.(x„E) 



det 



'd^s,{^,Ey 



-1/2 



(29) 



where 



D(q+,q-,E) = (-l)^det 



92 s 


92 s 


9q~9q+ 




92 s 


92 S 




5£;2 



(30) 



The details of the derivation leading to Eq. (p9|) can be found in Appendix ^ 

For time-reversal invariant systems, to every j-trajectory on C going from q~ to q^ 
(solution of the integral J^^) there is a corresponding time reversed pair going from q+ to 
q^ (solution of the integral lj\ Evidently both terms contribute with the same stationary 
phase and amplitude to the sum in Eq. (p5|). Hence 



C,(q+,q-,E) 



2A 



(rf+i)/ 



— 5:e-^*^/^|D,(q+,q-,E)|^/2cos 



S,(q+,q-,E) 



+ 7i + ^0 



(31) 



where = sgn[5^S'j(x, E)/dp'^]. This is essentially the main finding of Ref. 
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Let us now discuss the conditions under which Eq. (^) fails. When q~ approaches q"*", 
the points xj, represented by astericks in Fig. ^, move closer to the energy shell. As we 
learned in Section ^ this is a case where the semiclassical approximation for the spectral 
Wigner function, from which Eq. (|3l[) was obtained, fails. Remarkably, even in this limit 
where q~ —>■ q^ there are "long" chords whose center points Xj, indicated by dots in Fig. |^, 
are typically far from C. Such contributions to Ce are still well described by Eq. (pTj), but 
are obviously unrelated to the "short" chords. The approximation scheme for W developed 
in the previous section to deal with the classical contributions due to the "short" chords 
centered at points x^- close to C, can be also used to obtain C^- In the remaining of this 
section we pursue this path for two very different energy smoothing regimes and obtain a 
semiclassical approximation of valid for any spatial separation scale within the classical 
allowed region. 

At this point an important remark is in order. The chord structure for a fixed distance 
|q+ — sketched in Fig. ^ is robust upon changes of q = (q"*" + q~)/2 unless one of the 
points q^ reaches the boundary of the classical allowed region. When this condition is met, 
the "short" chords coalesce with the "long" ones, corresponding in Eq. (^) to the case of 
coalescing stationary points. A semiclassical investigation for a similar situation was already 
reported [^], but it did not address wave function correlations. This is a most interesting 
physical situation because of its relation to tunneling rates and possible implications to the 
already mentioned Coulomb blockade systems []26|-pO[]. Unfortunately we were not able to 
develop a semiclassical approximation to this problem yet. We avoid it in this paper by 
restricting our analysis to points q^ inside the classical allowed region and distant from its 
boundary by a couple of de Broglie wavelengths. 



A. Large e ^ eiargc smoothing regime 

For e ^ ^largo the cutoff parameter e suppresses all but the shortest trajectory contri- 
bution to Eq. (0). The latter connects both tips of a "short" chord (see Fig. As 
discussed before, when the center xq of approaches the energy shell C, we have to employ 
the corresponding uniform approximation to W, given by Eq. (PP]). 

Instead of directly Fourier transforming the semiclassical spectral Wigner function, it is 
advantageous to step back, use W as given by Eq. (^) and invert the integration order. 
That is 



C,(q+,q-,i5;) ~ — - / rft e-^l*l/^Fo(q+,q-,t)exp 



iEt 
IT 



(32) 



where 



Fo(q+,q-,t) = 



{2TiTiy 



dp 



det[U + A^o(x,t)] 



5o(x,t) + p-(q+-q-) . (33) 



The above integral is evaluated by the stationary phase method. The stationary phase point 
Po = Po(q''', ^ ~ )i) is the solution of 



€oq — 



d_ 
Sp 



5n 



q+ + q- 



q 



(34) 
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In analogy to the case that lead to Eq. (]28|), the phase factor is best written in terms of the 
standard text book action R, with variables in the configuration space, namely 



R{)(q^, q , ^) = 'S'o(xo, t) + Po ■ (q^ - q" 



(35) 



Substituting the obtained stationary phase approximation of Fq into Eq. ( |32D we write Cs 
as 



dt g{t) exp 



where i'o{t) = sgn[9^S'o(x, t)/(9p^]. The function g{t) gives the amplitude 



(36) 



g{t) = 2'^ e-^l*'/^ 



det 



U + A^o(x,t) 



det 



-1/2 



(37) 



x=xo 



with xq 



(q, po), q = (q"*" + q )/2 and po = Po(q^) q , The phase $ stands for 

<l>(t,0)=Ro(q+,q-,t) + ^t, (38) 



where we introduced oc |q''' — q | as a control parameter of the integral, i.e. Xq ^ C when 

The integral in Eq. (|36D is dominated by the stationary phase points located at ±i(:o(6', -E) 
corresponding to the traversal time of the shortest trajectory going from q~ to q"^ and the 
one running backwards in time. As 0, ±to{0,E) coalesce at the origin. This situation 
again is very similar to the one encountered in the previous section when dealing with the 
spectral Wigner function for x near the energy shell [i.e. that one that leads to Eq. (pO|)]. 
The difference is in the functional form of the phase $ and on the behavior of g{t) near the 
origin. For a Hamiltonian of the form Eq. (^) we show in Appendix ^ that both the phase 
$ and the amplitude g{t) have a singularity at the origin. Notwithstanding, the integral in 
Eq. ( PB| ) is finite and can be evaluated using the uniform approximation method The 
result is 



C,(q+,q-,E) 



27rA 



{2TTh) 



+1 



\Do{q^,q-,E)\ So(q+,q-,i?) 



1/2 



Jd 



So(q+,q-,i?) 



(39) 



where we left out the smoothing factor e~^l*ol/^, since for practical purposes the condition 
to ^ ^/^ is always met. Details about the evaluation of the integral in Eq. (|36D leading to 
(^) are found in Appendix |B[ 

The above semiclassical approximation for Ce(q+, q^, E) is valid for any separation |q"^ — 
q~ I , provided the arguments belog to the classical allowed region. Let us examine the small 
and large separation limits. As q~ — > q"*" the shortest trajectory on C is well approximated 
by the "short" chord ^q, and the action So turns into 



So(q+, q , ^) ^ Po(q) |q^ - q" 



(40) 



where po(q) = j2m[E — V"(q)]. Hence, the determinant Dq simplifies to 
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|Do(q+,q",^)| 



1/2 



m 



Po(q) 



{d-3)/2 



q" 



\(d-l)/2 



(41) 



(The derivation of Eqs. 
write 

C(q+, q- 



I) and ( PI) is found in Appendix 0.) Collecting the results, we 

{2nY/^ m p^qY''^ -^d/2-1 [p(q)|q+ - q"IA 



i27rhY p^{E) 



p(q)|q+ - q^\/h 



d/2-1 



(42) 



and recover, by an appropriate normalization. Berry's original result Eq. (0). In the corre- 
sponding opposite hmit, when |q^ — q^| is large, we use the asymptotic expansion of the 
Bessel function for large arguments, namely, Ju{x) ~ yj2/{TTx) cos{x — 1/71/2 — 7r/4) and 
retrieve the semiclassical approximation given by Eq. (|3TD for the shortest trajectory. 



B. Small e < eiarge smoothing regime 

As the smoothing parameter e is shrinked, longer trajectories have to be taken into 
account. As a consequence, when e < ^large, the approximation scheme becomes subtler 
than the simplified one discussed in Section [111 A|. Before exploring this regime, it is useful 



to remind ourselves that the semiclassical contributions to the spectral Wigner function 
come from orbits connecting the tips of either "short" or "long" chords, as depicted in 
Fig. 0. We thus classify the classical trajectories into three categories, (a) Trajectories 
connecting tips of "long" chords. We already discussed this case at the beginning of Section 
IT|. Here, since the trajectories are isolated, we just use Eq. (|31|). Difficulties, if any, arise 



from the large number of trajectories entering the sum, as regulated by the e. (b) Short 
trajectories connecting "short chords" . Those were analyzed in the foregoing subsection, (c) 
Long trajectories connecting "short chords". The analysis of Ce(q^,q~,£') becomes much 
clearer in the two cases discussed in the following, namely, x close and far from the energy 
surface C. 

When the long trajectories have tips on "short" chords, and q~ approaches q"*", the 
centers of these chords approach C. Here in addition to the trajectories of type (a) and 

(b) , we have to account for the long trajectories associated with "short" chords of the type 

(c) but whose centers are close to C. Hence, as already discussed, for evaluating points close 
to the energy shell we must use a suitable W for Eq. (^. This is just Berry's scar expansion 



formula in terms of periodic orbits [11,32,33|. Hence, for small distances |q^ — q |, we write 



c,(q+, cr,E)= Yl c'^lq^^ q". ^) + ^) + (q+, ct.e), (43) 



jg"long" 
chords 



where C| stands for the j-th type (a) orbit corresponding to a term in the sum in Eq. (plD, 
C° for the type (b) trajectory term given by Eq. (|55|), and C^^^^ is the result of Eq. (^ 
for W taken as the scar expansion formula. We do not intend to present here any detailed 
analysis of this integration, but it is easy to realize that the final result will be an expansion 
in terms of all the periodic orbits that pass through the points q^. Such conjecture is 
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further supported by noting that when q = q"*" = q the integration Eq. (|D reduces to the 
projection 'down' p of the Berry's scar expansion. The latter case corresponds to Bogomolny 



formula I^S] for scars in the probability density in configuration space, which involves all 
the periodic orbits that pass through the point q (see Ref. |]3^). 

In the other limit, when |q^ — q^| is large, the centers Xj of all the chords sketched in 
Fig. 1^ are far from C, and we thus write 

C,(q+, cr,E) = Y: C^(q^, cr,E) + C°(q+, q", , (44) 

where we separate the contribution of the shortest trajectory given by the uniform 
expression Eq. (p^). 



It is interesting to note that both formulas Eq. (^) and (^i]) imply that the spectral 
autocorrelation function has contributions arising from W taken at points well inside the 
energy shell. In the case of Eq. (|43| ) these are the trajectories with tips on "long" chords. In 
general, these contributions can not be neglected, as already discussed at the end of Section 
|I[ Particularly, the "off-shell" scars of the spectral Wigner function provide the only 
periodic orbit contributions for in the case of large separation |q+ — q~|. 



C. Spatial averaging with e < eiarge 

We now investigate the effect of spatial averaging on the spectral autocorrelation func- 
tion. We are interested to know under which circumstances this averaging washes out the 
system specific features, allowing one to define a local system independent (or universal) 
regime for Cg. Furthermore, the additional spatial averaging brings our results into close 
relation to numerical experiments, such as the ones in Refs. [pi|,p^. 

Defining the local spatial average as 

(C,(q+,q-,E))q= JjciCe{ct,^-,E) , (45) 

where f2 is a configuration space region of surface A{Q) (in d-dimension) , covering many 
"local" de Broglie wavelengths Aq across. We define Aq = 2'n'h/p{q) in terms of the local 
momentum p(q), consistent with the semiclassical approximation. The average is restricted 
to a region Q of classically small variations of the smooth potential ^(q). 

The ubiquitous robustness of Berry's expression for the autocorrelation of chaotic wave 
functions can be attributed to the following: At small separations |q^ — q^| the local space 
average kills the terms C| associated with the "long" chords. This suppression is due to 
the fact that such terms oscillate in a scale smaller (or at most comparable) to the de 
Broglie wavelength. Indeed, for small separation the "long" chords are almost parallel to 
the momentum space and thus its centers have pj(q+,q^,i?) ^ (see Fig. /^From 
Eq. (p8|) , Sj(q^, q^, -E) ~ Sj{yij,E) and thus the local spatial wavelength in Eq. ( ^Tf) is 
approximately 



Ajq 2-Kh 



dSj{^j,E) 



1^ . («) 

I s jp I 
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This is just about the de Broghe wavelength since \^jp\ ~ 2p(q). Moreover, the requirement 
of small spatial averaging regions assures that p(q) is approximately constant for q inside Q. 
Consequently, C° for small Iq"*" — q~ | coincides with Berry's result and remains unaffected by 
the local spatial average [^. Furthermore, if none of the wave functions inside the energy 
window from which Ce is built shows a strong visual scar due to periodic orbits, which is 
often the case, there is no reason to expect a sizeable correction due to {C^'^^^)^- Thus, 
after the averaging the leading contribution for Cs would be the local system independent 
expression Eq. (^). 

Extrapolating our semiclassical analysis of (Ce(q^, q^, -E'))q to energy smoothings of the 
order of one energy spacing {i.e., e ~ ^smaii), we find that our results are still consistent 
with the numerical investigations |plip2| of the autocorrelation function C(q+,q^) [in this 

- q~|. For instance, in 
are small unless the 



case Eq. (0)] on individual eigenf unctions for small values of |q 
Ref. ||2^ we observe that the corrections to Eq. (^) given by {Cl 
eigenfunction has a strong visual scar due to a simple periodic orbit. Likewise, as Eq. 



symmetric with respect to the orientations of (q^ 
observed angular dependence of C(q^ 



— q ) for a fixed value of 
arrives precisely from {Cl 



q 



§ IS 

the 



,q~) arrives precisely from {C^^^'^)^. 

Before concluding we like to add that the local spatial average is the mechanism respon- 
sible for the elimination of the contributions of trajectories associated with "long" chords in 
the Bogomolny scar formula for the spatial probability density |35|. In our formalism, the 
latter is recovered after making the local average Eq. (|45|) over the autocorrelation function 
Cs, Eq. (|43|), and then taking q+ q^. In other words, the Bogomolny scar formula cap- 
tures (in configuration space) only the scar contributions of periodic orbits near the energy 
shell of the spectral Wigner function, since the "off-shell" scar contributions are washed 
out by the local spatial average. 



IV. CONCLUSIONS 

We investigated the spatial two-point autocorrelation of energy eigenfunctions '?/'ri(q) 
corresponding to classically chaotic systems in the semiclassical regime. We use the Weyl- 
Wigner formalism to obtain the spectral average Ce(q"'", q~, E) of '?/'„(q"'")'?/'*(q~), defined as 
the average over eigenstates within an energy window e centered at E. In the considered 
framework Ce is just the Fourier transform in momentum space of the spectral Wigner 
function iy(x, E; e). 

The advantage of this formalism comes from the observation that W is almost like 
tailormade for semiclassical approximations. At each phase space point, x = (q, p), the 
semiclassical behavior of W is associated with all the classical trajectories on the energy 
shell E whose end points are joined by a chord centered at x. These classical contributions 
are exponentially suppressed when the trajectory traversal time is t ~ h/e. 

In distinction to most studies so far, this paper addresses smooth Hamiltonian systems 
of the form given by Eq. (^. Our results can evidently be straightforwardly employed to 
calculate Ce in billiard systems, where the phase space structure is much simpler than the 
one considered here. For smooth systems, we show that it is still possible to distinguish 
in Ce between a local system independent regime and another one that carries the system 
classical chord structure information. This is obviously also valid for billiards. The interplay 
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between the spectral average window, which controls the upper time scale of the classical 
contributions, and the spatial separation scale dictates which aspects prevails. As a result 
we obtain semiclassical expressions that bridge the existing formulae for the autocorrelation 
function C^. 

To the best of our knowledge, the studies of Cg found in the literature are based on 
Green's function methods, and employ a given arbitrary separation between "zero-length" 
and "long" trajectories. In billiards, due to their simplicity, the semiclassical "zero-length" 
Green's function is an excellent approximation to calculate C^, even for e comparable with 
eiarge, provided that the separation is not too large. In smooth systems corrections account- 
ing for the energy surface curvature become rapidly necessary as the spatial separation is 
increased. Such corrections, albeit in principle feasible to obtain with the Green's function 
method, are easier to estimate with the here employed framework. This is an important 
advantage of the formalism we employ. Our study goes beyond that issue showing that 
the Wigner-Weyl formalism is a quite general framework for semiclassical approximations, 
clearly revealing the inextricable relation between the classical chord structure and the 
choatic wave function correlations. 
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APPENDIX A: DEMONSTRATION OF EQ. (28) 



In this Appendix we show that both amplitudes in Eq. (^) are identical, thus proving 
that Ref. addresses one of the limiting cases of Cg studied in this paper. For the sake 
of clarity, it is convenient to express both sides of Eq. (pOf ) in terms of the action defined 
as R(q+,q^,t) = S(q+,q^,£') — Et. The variation of R with respect to its independent 
variables, namely, q''',q~, and t gives 

= ±dR/dci^ and -E = dR/dt . (Al) 

It is also convenient to introduce a short notation for the second derivatives of R 

_ a^R ap+ ^ _ a^R ap- 
R4-+ = — — -7- — — — - Rj = 



9q+^ 9q+ (9q+9q~ 9q+ 

R = R ^^^J_^^ .A2) 

^ 9q~5q+ 9q~ (9q~^ dc\^ 

which form a set of four (i-dimensional matrices. 

With the elements at hand one readily writes D{c\^ , q~ ,E), as defined in Eq. ( pOD by a 
{d + 1) X [d + l)-dimensional matrix determinant (see, for instance, Ref. ^) 
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D = (0) ' det(-R^^) = - det(-R^_) . (A3) 

We switch now to the Lh.s. of Eq. (PUj). The 2d x 2(i-dimensional stabihty matrix A1(x) 
is related to the second derivatives of the action R as follows 

— R_i_^R — -^+- 

^R_| R — R-| ^R_| 



^(x) = . (A4) 



The justification for the above equation can be found in the Appendix A of Ref. [|ri|]. It is 
now straightforward to show that 

1 , . y/ M detfRj — — R++ — R — + R — u) , . ^, 

dat[a+A<(x)]^ det(-R,_) '^^^ 

It remains only to express d^S/dp'^ as a function of R to conclude the demonstration. 
The relation between the second derivatives of the center action ^(q, p, E) and the stability 
matrix A4{x.) can be obtained by differentiating q"*" = q'''(q^,p^) and p"*" = p'''(q~,p~) 
with respect to q and p. By writing q~ = q~(q, p), p^ = p^(q, p), and recalling that 
q^ = q =F {dS/dp)/2 and p^ = p ± {dS/dq)/2, we arrive at 



[U - JBi^)] = A^(x) [U + JBi^)] , (A6) 



where 



= 2 ii J=[\ I) ■ (AT) 

Using the obvious identity J'~^ = — J we recast -B(x) as 

S(x) = -M(x)][a + A^(x)]-^ (A8) 
The above relation with the aid of Eq. ( [A^ ) renders 

= 22'^(R+„ - R++ - R__ + R_+)-i . (A9) 

Finally, by collecting (|A3|) , (|A5|) , and ( [A9| ) we arrive at 

-(§){det[. + A.(x)]det(0)p^ 
_ _ ldt\ det(-R^_) _ D_ 

- "22'^' ^^^^^ 

which proves Eq. (^). 
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APPENDIX B: DERIVATION OF THE UNIFORM APPROXIMATION EQ. (M) 



In this appendix we evaluate the integral Eq. ( |5B| ) by the uniform approximation method, 
to obtain Eq. (pQ]) for the spectral autocorrelation function. As already mentioned, the 
integral in Eq. (^) is dominated by its stationary phase points ±tQ{9,E) (solutions of 
dKo/dt + E = 0), that coalesce at the origin as ^ 0. Thus, we start analyzing the 
structure of the phase 9) and the amplitude g(t) near t = 0. 

For 6) we replace the action Rq by the one found in Eq. (|35D. We then expand the 
center action 5'o(x, t) [providing that po = Po(ct^5C[ is fixed by the condition Eq. (0)], 
up to third order in t, as in Eq. (JT^). If the system Hamiltonian has the form of Eq. 
the expansion reads 



so for Eq. (|34D we have 
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dV 


' 1 d^V 


fi 

\ m 


(9q 


+ 2? A 2 



p 



(Bl) 



dp 



-A{q,t) Po = (q^ 

m 



q 



(B2) 



where Xo = (q, po), with q = (q"*" + q )/2 and A is the d x (i-dimensional matrix 



A{q,t) = 11 + 



■(q) 



Solving Eq. 



12m (9q^ 

for Po and using that A^^ = U — ■^^d'^V / dc^ + 0{t^), we arrive at 

t d^V 



(B3) 



Po(q ,q ,t) 



m 
T 



iq -q 



12 dcf 



(q) + Oit' 



(B4) 



where m(q^ — q )/t comes from the first order term in the t-expansion of the center action 
Eq. (pTl). Eqs. M and (M§ yield 



$(t,^^) = ||q^ 



q 



[E-Vici)]t + ^t, 



(B5) 



Here the explicit part of $ arises from the lowest order in Eq. ( plD and <l>(t, 6) is an analytical 
function of t. The phase ^(t, 6) has a pole of order one at the origin of the complex plane t. 
The behavior of g{t), defined by Eq. (^), near t = is dominated by 



^1/2 



m 



d/2 



x=xo 



U\d/2 



det 



^(q,t) 



-1/2 



(B6) 



as seen from Eq. ( p^) . For clarity purposes we introduce g{t) an analytical function of t, 
defined as g{t) = g{t)/\t\'^/^. 

The center action 5*0 and Po(q^,q~,^) are odd function of t ||33|, and consequently so is 
the phase 9). Since g{t) is an even function of t, we are allowed to express Eq. (|36|) by 



18 



2A 



(27r^)^/2 



- Re I{h, 9) 



2A 



2+1 



Re 







dt ^iri- exp 



(B7) 



where we used that I'oit) = sgn[c?^S'o(x, t)/9p^], which is the difference between the number 
of positive and negative eigenvalues of 9^S'o(x, )f:)/(9p^, is simply equal to —td/\t\ when the 
energy shell is convex (the case considered in this study). We have now to deal with a single 
stationary phase point to{9,E), that coalesces with the lowest limit of integration, and for 
which the stationary phase reads 

$(to,^) = Ro(q+,q-,to) + ^to = So(q+,q-,^) • (B8) 

The uniform approximation to Eq. ( [B7| ) involves the suitable change of the integra- 
tion variable, w = w(t) : [0, +oo) — > [0, +oo), such that it is invertible and reduces the phase 
of the integrand to the canonical form 



1 

2 



w 



w 



(B9) 



adhering to ^{t,6) as given by Eq. ( P5|) . Likewise, we must require that the stationary 
points ±to of $ correspond to the saddle points ±wo = ±z{9) of {i.e. dcj) / dw\w=±ujQ = 



for ±wo = w{±to)). This is achieved by making z{6) = So(q''',q 
after the change of variable the integral in Eq. (p7|) becomes 



iih,9) 



^ G{w) 
aw — — exp 



w 



d/2 



where 



G{w) m 



dt 



id— 
4. 



\w 



d/2 



|i|d/2 



dw 



E) [see Eq. (|B8|)]. Thus, 
(BIO) 

(Bll) 



and G is an even analytical function oiw. In general the next step of the method of uniform 
approximation is to expand G around wq [^, but in our case is sufficient to keep only the 
first term, namely 

dt 

G{w) ^ G{wo) = g(to = t{wo)) — \wo 



where 



dt 



dw 



\d(j)/dw\ 



10=100 



dw 



1 



d/2 



W = Wo 



|So(q+,q-,i?)|V2 



dto 



dE 



1/2 



\d^/dt\u=,. 

The last equality was obtained by applying the L' Hospital rule. Now, recalling Eq. 
write 

G{wo) = e-^l*«l/^|Z}o(q+, q", E)\'/' So(q+, q", EY'^-'y' , 
and using an integral representation of the Hankel's function H^\z) we have 



(B12) 



(B13) 



we 



(B14) 



vr G{wq) 



H 



(1) 

d/2-l 



So(q+,q-,E) 

h 



So(q+,q-,i?)'^/2-i 

Finally, collecting the results of Eqs. ( pi5|) and ( pi4|) and taking the real part in Eq. 
we obtain the uniform approximation to Eq. (pOf). 



(B15) 
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APPENDIX C: DEMONSTRATION OF EQS. (40) AND (ET 



In this appendix we derive Eqs. ( PP] ) and (PT|), approximate expressions for the action 
So(ci^) q~) and the determinant Dq respectively. Both relations are valid provided the 
control parameter 6 oc |q^ — Q I is small. Eqs. (|40|) and (^) show that our result for the 
spectral autocorrelation function Eq. (|39|) reduces to Berry's one. 

Let us start with Eq. (|40|) . The action Sq is given by Eq. (B8) where is determined by 
the stationary phase condition dKo/dt = —E. For 6 small, Rq is 



Ro(q+,q ,t) 



m 
2t 



|q -q 



r(q)t 



(see Eq. (|B5| ) and the discussion preceding it). The stationary phase point to is thus 

m . , 



tn 



|q -q 



(Cl) 



(C2) 



where Po(q) = y2m[E — V^(q)]. By using Eq. ( |Cl|) and Eq. (|C2|) in Eq. ( pSj) we arrive to 
the approximation Eq. (^Ol). 

For the determinant Dq, from Eq. ( [AlOp we write 



\Do{ci^,q-,E)\ 



1/2 



dto 


1/2 


dE 





det [U + A^o(x,to)] 



-1/2 




-1/2 









(C3) 



where xq = (q, po) is the center of the shortest chord with tips on a classical trajectory 
(see Section [III A]) . We used the center action S'(x, t) instead of S'(x, E) because for any of 
both ^q(x) = —dS/dp. When 6* — 0, the center xq approaches the energy shell, the chord 
£o ~ toi^ [to is given by Eq. (|C^ )] well approximate the classical trajectory in phase space, 
and hence becomes the identity map, yielding 



dto 


1/2 


d'Ro 


-1/2 


dE 




dt^ 





m 



t=to 



Po(q)^/^ 



|q - q 



-1 1/2 



(C4) 



where we used Eq. ( |C1|) for Rq. Up to the same order considered in Eq. ( |Cl|) the matrix A 
in Eq. (|B6|) becomes the identity and 



det 



' d^So{^,to 



-1/2 



m 



d/2 



Po(q)'/' 



x=xo 



|q^ - q 



\d/2 



(C5) 



Collecting these results in Eq. (|C3) we arrive at Eq. (HI] 
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FIGURES 




FIG. 1. Chord structure of the semiclassical Weyl propagator, Eq. (0), and of the spectral 
Wigner function, Eq. (13), in a 2(i-dimensional phase space. For a central point x we show typical 



chords with their tips connected by classical trajectories on different energy shells but with the 
same traversal time t. The dashed line indicates the locus, on a given energy shell E, of all the 
tips of chords centered at x. The classical trajectories connecting these type of chords are the 
semiclassical contributions to W. The dashed area represents the symplectic area corresponding 
to the action 52(x, i?) (see text). 




FIG. 2. Illustration of a composition of trajectory segments that closely follow a periodic orbit. 
Berry's scar formula single out the trajectories connecting the tips of the two depicted chords which 
become indistinguishable from the period orbit as x approaches C. 
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FIG. 3. Chord structure of the semiclassical approximation for the spectral autocorrelation 
function C^. The horizontal dashed lines represent the projection on configuration space of the 
chords matching the vectors ±(q"'" — q~). The points of stationary phase are located on the 
d-dimensional associated momentum space to q = (q+ + q^)/2 and indicated by astericks (*) and 
dots (•) (See text for details). C stands for the (2d — l)-dimensional surface of constant energy 
E, whereas the axis q and p represent d-dimensional surfaces. The semiclassical contributions to 
Ce are given by all trajectories on C flowing between the tips of chords ^q. Panel (a) illustrates a 
typical time reversal symmetric case, while (b) a case when this symmetry is absent. 



24 



